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Abstract. A fundamental requirement in the new generation of high resolution Cosmic Microwave Background 
imaging experiments is a strict control of systematic errors that must be kept at pK level in the final maps. 
Some of these errors are of celestial origin, while others will be generated by periodic fluctuations of the satellite 
environment. These environment instabilities will cause fluctuations in the measured signal output thus generating 
correlated effects in the reconstructed maps. In this paper we present an analytical study of the impact of periodic 
signal fluctuations on the measured sky maps produced by the Planck survey. In particular we show how it is 
possible to estimate analytically the damping factor of the peak-to-peak amplitude of the fluctuation at the 
instrument output after the projection in the final maps. 
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, ^ 1. Introduction 



The anisotropy pattern in the Cosmic Microwave 
Background (CMB) radiation carries information about 
fundamental cosmological properties such as the present 
expansion rate and average density of the universe, the 
amount of dark matter, and the nature of the seed fluc- 
tuations from which all structures in the universe arose 
(Vittorio & Silk |1984 Efstathiou |1987D . In 1992, the 
COBE-DMR maps provided the first ex perim ental evi- 
dence of such fluctuations (Smoot et al. 1992 ) at angu- 
lar scales > 7° at a level ST /T ^ 10"^ (where T = 
2.725 ± 0.002 K, Mather et al. |1999D . Since then, several 
balloon-borne and ground-based experiments have pro- 
vided CMB anisotropy detection at sub-degree angular 
scales, thus improving substantially our knowledge of the 
anisotropy power spectrum. Recently, sub-orbital experi- 
ments have been able to map small sky patches with high 
angular resolution and sensitivity (e.g. Netterfield et al. 
20001, De Bernardis et al. |200C|, Lee et al. pOOll Stompor 



where the entire sky is in full view and the environmen- 
tal conditions are much more favourable (Danese et al. 



1996). The NASA satellite MAP (launched in June 2001) 



et al. pOOl| , Halverson et al. |200l[ Pryke et al. g00l|). 

Full-sky CMB maps with high resolution and sensitiv- 
ity, however, can be obtained only by space experiments, 
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will provide the first improved full-sky maps since COBE, 
while the ESA satellite Planck will map the CMB with 
an accuracy set mainly by astrophysical limits. 

A key factor for the success of these missions is a 
strict control of systematic effects (see, e.g., Bersanelli & 
Mandolesi, 2000), particularly those that are synchronous 
with the spacecraft spin period. These fluctuations, in 
fact, may be particularly difficult to remove from the data 
stream because of their similarity with the true sky signal, 
although some data analysis methods have been developed 
to deal also with spin synchronous signals of celestial ori- 
gin (Delabrouille et al. ^00l[ ). 

In many cases the above effects display a periodic be- 
haviour leading to a periodic spurious contribution to the 
measured signal. While the transfer function from the fluc- 
tuation source to the measured signal depends on the de- 
tails of the instrument and of the effect being considered, 
the transfer from the measurement output to the sky map 
does not depend on the source but only on the period of 
the effect and on the scan strategy, and therefore it can 
be treated independently. 
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This work is a study of the impact of periodic signal 
fluctuations on the measured sky maps to be produced 
by the Planck survey. In particular we show how a pe- 
riodic fluctuation in the instrument output transfers to 
the flnal maps and analyse the behaviour of the damping 
factor versus the oscillation frequency (Sect. In this 
section we also evaluate the possibility to reduce the level 
of residual systematic errors by applying destriping algo- 
rithms (Sect. 2.2.2), also considering the combination of 
such systematic errors with instrumental noise. The main 
result of this analysis is summarised by a general relation- 
ship that allows us to estimate the final peak-to-peak error 
in a map starting from a periodic effect of arbitrary shape. 
In Sect. ||we discuss the application of this relationship to 
an example of particular interest for Planck-LFI. Finally, 
in Sect. ^ we draw our main conclusions. 

Although the examples presented in this paper are rel- 
ative to Planck, the concepts at the basis of our treat- 
ment are general, and can be applied to any CMB imaging 
experiment. 



2. Transfer of periodic oscillations from 
instrument output to maps 

2.1. Level of residual fluctuations in scan circle 
time-scale 

The Planck measurement strategy is such that the signal 
coming from each resolution element in the sky will be 
measured approximately 60 times during one scan, thus 
reducing the effect of fluctuations occurring on time-scales 
greater than the spin period. 

In particular, if we average N measurements of the 
same pixel taken at times tj then the residual signal os- 
cillation is given by: 



1 ^ 

(<5T,ky(i))t„ =-^<5rsky(t,) 



(1) 



data stream after averaging over N scan circles is simply 
given by: 



Af sin(7riVTspi„/Tf ) 



N sin(7rTspin/Tf) 



(3) 



Two special cases can be derived from the above equa- 
tion: the case of "spin synchronous" fluctuations (for 
which Tf = Tgpiii/fc where k is any integer) and the case 
of "spin resonant" fluctuations (for which T{ = N Tsp-„Jk, 
where k is any integer not multiple of N). 

Spin synchronous fluctuations are not damped by the 
measurement redundancy^ while spin resonant fluctua- 
tions, instead, are such that for every pixel the average 
fluctuation after N consecutive measurements is zero. 

Let us now consider first oscillations with Tf < rgpin- 
In Fig. |l| we show a plot of {5T^-^^)n given by Eq. ^ for 
values of Tf up to 70 s, considering Planck nominal values 
(^spin = 60 s and N — 60) and assuming Af — 0.5 mK. 
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Fig. 1. Plot of ((5rP^yP)jv for for Af = 0.5 mK and Tf 
between 20 s and 70 s. The graph shows the two peaks 
relative to the spin synchronous periods of 60 s and 30 s 
(for which {5T^]^^)n = 2Af) and a series of spin resonant 
periods (for which {5T^\^^) n = 0, see inset). 



where tj — t + {j — 1) Tgpin (Tgpin is the spin period which 
is also the interval between two consecutive measurements 
of the same pixel in the sky) . 

Let us consider a sinusoidal signal oscillation with am- 
plitude Af and period Tf, i.e. 5Tsky(i) = v4f cos(27rt/Tf). 
Eq. (|l|) takes then the form: 



(<^r,ky(<))w = ^^cos(27rt,/Tf) 



N 



\ 1^, , (AT 1^ \/ -I sm(7riVTspi„/Tf) 

= cos[7r(2t+(Ar-l)r i„)/Tf]^^ 

N sm(7rTspi„/Tf) 

Eq. (H) gives the amplitude of the signal oscillation 
after averaging N measurements of the same pixel starting 
at time t. Because the term cos [7r(2t -I- (A^ — l)Tspin)/Tf] is 
bounded by ±1, then the peak-to-peak fluctuation in the 



The picture shows that the residual peak-to-peak 
fluctuation in a sky scan circle equals 2Ai for spin- 
synchronous periods and is of the order of Af/N other- 
wise. The inset in Fig. ^ highlights that for spin resonant 
periods the oscillation is completely cancelled by the av- 
eraging. 

Fig. ^ shows the case of oscillations with Tf > Tspin- 
In this case there are no more spin synchronous periods, 
but we still have an effect of resonant periods up to Tf = 
N X Tgpin (which equals 3600 s in our model) . The envelope 
defining the upper hmit of {ST^i^^) n (solid fine) is given 
by: 



/N 



Ai 



1 



N sin(7rTspin/Tf) 



(4) 



^ This implication can be easily derived from Eq. consid- 
ering that lim^;-,! (sin(A'"7ra;)/ sin(7ra;)) = N 
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Fig. 2. Plot of {ST^i^^)N for for Af = 0.5 mK and n 
between 60 s and 10000 s. The graph shows that no more 
spin synchronous periods exist after rt — 60 s, while the 
last spin resonant period is at Tf = 3600 s. The solid line 
represents the envelope defined by Eq. (0). 



Some remarks need to be made concerning the rele- 
vance of spin-resonant and spin-synchronous periods in 
the context of Planck measurements. 

Spin-resonant fluctuations in principle do not con- 
tribute to the overall effect; in reality many non- ideal be- 
haviours (like instability in the spin period, pointing er- 
rors, variations in the repointing time, etc.) will concur 
in eliminating such sharp resonances. The impact of spin- 
synchronous periods (in particular T{ = 60 s) , instead, will 
be less sensitive to smearing by scan non idealities. In fact 
if we consider the main peak in Fig. |l| we see that it has a 
width (at half maximum) of ~ 1 s. Because the Planck 
spin rate will be stable at the level of ^ 10~^ rpm/h it 
follows that even in the case of a constant drift in the spin 
rate a fluctuation which is initially spin-synchronous will 
maintain this characteristic for time-scales of the order of 
one week. 

Therefore, in our paper we will adopt a conservative 
approach by neglecting only the effect of spin-resonant 
periods. This implies that Eq. (|^) is used to estimate the 
effect on the measured data stream of fluctuations with 

"H" ^ '^spin- 

2.2. Level of residual fluctuations in final maps 

The simplest way to estimate the impact of periodic fluc- 
tuations on the final measurements is to project the data 
stream of measurements averaged over multiple scans onto 
a map with a given pixel size and then evaluate the peak- 
to-peak amplitude. This rather simple procedure gives a 
gross upper limit of the amplitude of the effect, which is 
usually greatly reduced by applying so-called "destriping" 
algorithms to the time ordered data. 

In the following we derive a general relationship to 
estimate the final peak-to-peak amplitude of a given sys- 
tematic effect by knowing its amplitude at the instrument 
output, its frequency and the map pixel size. This rela- 



tionship is obtained analysing the simple projection of the 
periodic signal onto the sky map and then considering the 
additional damping factor obtained after the application 
of our destriping code to the data stream. 

2.2.1. Projection of averaged data streams onto sky 
maps 

The projection of a time ordered array of data onto a map 
consists in the association of every data sample to a pixel 
coordinate in the sky. 

If we consider the Planck scanning strategy it is clear 
that pixels close to the ecliptic plane will be "visited" 
by a smaller number of scan circles compared to pixels 
around the ecliptic poles. Therefore, the amplitude of a 
fluctuation projected onto a sky map will be maximum 
along the ecliptic, where each pixel is crossed by iVscan ~ 
6'pixei/^rep circlcs, whcrc 0pixoi is the pixel size and 6'i.op 
indicates the repointing angle, that in Planck current 
baseline is ~ 2.5'. 

A rather simple analysis of the data averaging in pix- 
els around the ecliptic plane leads to the following formula 
for the peak-to-peak amplitude of a fluctuation with am- 
plitude A[ and period rt when projected onto a sky map 
with pixel size ^pixci: 



W7) 



At 



sky /map 



N X Ns, 



sin(7riV X A^'scanTspin/rf) 



sin(7rTspin/Tf) 



(5) 



From Eq. ^ it appears that {ST^^^^)n^ap can be zero 
for every value of Tf that satisfies the relationship ti — Nx 
-^scan'^spin/fc, with k integer. A set of these values is what 
we have called spin resonant periods, i.e. those periods 
which are resonant with the scan circle time N x 



spin • 



Another set of these values corresponds to periods that 
are resonant with the repointing time Tscan — -^scan x Tspin- 
Considering that A'scan is not constant all over the map 
(pixels close to the ecliptic poles are crossed by a greater 
number of scans compared to pixels close to the ecliptic 
plane) this latter set of periods does not actually produce 
a global zero-effect, but only in a limited number of pixels. 

Therefore, with the conservative assumption that spin 
resonant periods will not be relevant in practice, we can 
conclude that the best estimate of the peak-to-peak ampli- 
tude {ST^^"^ ) i-n&p can be well approximated by the upper 
limit of Eq. (^), i.e.: 



= 2 



'map 



\ max 
/map 



NxNs, 



1 



sin(7rTspin/rf) 



2At 



map 



(6) 



where -Fmap is the damping factor of the oscillation from 
the instrument output to the sky map. 

We can now compare the damping factors calculated 
by Eq. (||) with those obtained by projecting periodic sig- 
nals with different values of Tf onto sky maps. The maps 
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(produced according to the HEALPix hierarchical struc- 
ture, Gorski et al. 1998) are relative to the position of a 30 
GHz LFI radiometer and have a pixel size of 13.7'. All the 
maps have been produced assuming 10200 hours observa- 
tion and the Planck baseline scanning strategy, which 
implies an average integration time per pixel of ^ 46.7 s. 

Fig. 1^ shows the result of this comparison and high- 
lights a remarkable agreement between the damping fac- 
tors obtained directly from maps and calculated from 
Eq. (I). 
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Fig. 3. Comparison of the damping factor as a function 
of Tf calculated by Eq. (^) (solid line) and obtained from 
maps (black dots) (0pixei = 13.7', 6'scan — 2.5'). 
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Fig. 4. Application of the destriping code to a "pure" pe- 
riodic systematic effect. The upper map is relative to a 
sinusoidal fluctuation with a rf = 4000 s and A{ = 5 mK. 
After post-processing the time ordered data with our de- 
striping code the peak-to-peak effect on the map is re- 
duced by a factor ^ 34. 



2.2.2. Removal of periodic systematic effects from 
Time Ordered Data 

In this section we evaluate the residual systematic error 
on the maps after the application of a "destriping" code 
developed by the LFI team for reducing the impact of 1// 
noise fluctuations. 

This code generally uses averaged one-hour scan cir- 
cles (although it has the capability to work with the full, 
unaveraged data streams); because averaging acts like a 
low-pass filter, only the very low frequencies remain in the 
map as a residual effect. The basic hypothesis is that for 
each circle the residual effect can be well approximated by 
a constant additive level related to the mean level of the 
spurious fluctuation during the scan. The code optimises 
these levels by minimising the map temperature values in 
pixels observed by more than one scan circle during the 
mission. Further details of the co de implementation can 
be found in Burigana et al, 1997t| . 

In Fig. I we show an example of the application of 
this code to a pure periodic fluctuation with a period of 
4000 s and an amplitude (at the instrument output) of 
1 mK peak-to-peak. A comparison between the two maps 
shows that the code is able to reduce the impact of such 
fluctuation by a factor of ^ 34. 



So far we have considered the application of the de- 
striping code to time ordered data containing only the 
periodic fluctuation, i.e. with no noise present. We have 
also shown that the presence of white and 1// noise in the 
data stream does not change the ability of the code to re- 
move periodic fluctuations. To verify this we have applied 
our code to a data stream containing noise (white and 
1//) plus a periodic signal and to a second data stream 
containing only white and 1// noise. 

The noise is typical of LFI 30 GHz radiometers, 
i.e. white noise and 1// noise, with a spectral den- 
sity of 230 fiK-Hz^^^^ and a knee frequency of 0.1 Hz 
(rather conservative estimates of the expected perfor- 
mances 30 GHz LFI radiometers)]^. The periodic fluctu- 
ation is the same as in Fig. ^ 

By taking the difference between the map containing 
noise plus the periodic fluctuation and the map containing 
only noise (both after destriping) we obtain a map that is 
virtually indistinguishable (at the level of 10~^ /iK) from 
the bottom map of Fig. ^, which was obtained by applying 
the destriping algorithm to the periodic oscillation with- 



^ The reported spectral density is relative to the noise re- 
quirement of the 30 GHz LFI radiometers and it has been 
chosen as a conservative case. The goal value is of about 
180 /iK • 'Rz~^^^ with an expected knee frequency of 0.05 Hz 
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out noise. This indicates that the algorithm was able to 
recognise and reduce the effect of the periodic fluctuation 
even in presence of noise stream with a much larger am- 
plitude level. 

Our next step has been to evaluate the destriping 
damping factor for periodic fluctuations as a function of 
Tf . Fig. P shows how this damping factor increases approx- 
imately linearly with Tf for periods greater than the spin 
period^. It is important to underline that the exact form 
of the linear interpolation shown in Fig. |^ is dependent on 
the conditions under which the maps have been calculated, 
which are, in the case discussed here, a sampling rate of 1 
sample every ~ 12' and a map pixel size of 13.7' (typical 
conditions of LFI 30 GHz channels). 
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Fig. 6. Total damping factor from instrument output to 
final map for the peak-to-peak amplitude of periodic os- 
cillations having a period Tf > Tspin considering the de- 
striping damping as in Fig. |[ 

3. Application to a periodic fluctuation of 
arbitrary shape 

In this section we extend our analysis to periodic fluctu- 
ations of general shape and show some examples of its 
application. Under very general assumptions we can write 
the oscillation in Fourier series as: 



+ 00 

E 

j=-oo 



Aj exp{i2T:i>jt) 



(8) 



Fig. 5. Destriping damping factor for periodic oscillations 
versus Tf. Dots indicate the values obtained from maps 
while the solid line represents the linear interpolation 
-Fdcstr(Tf) = CTf + 6 with a = 0.1067 s'^ and b = 1.7992. 

If we indicate with -Fdcstr('''f ) = o-Tf + b the destriping 
damping factor, we can write the following general form 
for the total damping factor of the peak-to-peak amplitude 
of periodic oscillations from instrument output to final 
map: 



F{t[) > Finapl-Tf) X ^dcstr(Tf) 



(7) 



where i^map can be calculated by Eq. (^) . Note that Eq. 
provides a lower limit of the total damping factor rather 
than the exact value. This is a consequence of our as- 
sumptions regarding the beam position and the scanning 
strategy which are marginally conservative with respect 
to the destriping efficiency. 

In Fig. ^ we show a plot of the function i^(Tf ) consid- 
ering the destriping damping factor relative to the case 
shown in Fig. ||. 



For spin-synchronous fluctuations no damping of the effect 
is obtained from the application of destriping codes 



where Vj represent the harmonic frequencies. We now ap- 
ply the same procedure discussed in the previous sections 
to each term of the series in Eq. (||) in order to calculate 
the peak-to-peak amplitude after one-hour averaging and 
after the projection of the averaged time data stream onto 
a map. Considering that for each harmonic term we have: 



N 



^ exp(i27ri/j(io + {k - l)Tspin)) = 



k=l 

= A, exp [ZTT (2io + {N~ l)r,^,^) v,] (g) 

Sm^TTTspini^j) 

then the final peak-to-peak amplitude for a general signal 
on the map after destriping can be written as: 



.ky /map ^ ^ ^ 



+ 00 

E 

j=-oo 



Aj / Fticst-c{Vj) 



sin(7rTspi„j/j) 



(10) 



Now we show an example of the application of this 
formalism to the signal oscillation shown in Fig. ^ both 
in time domain (upper graph) and Fourier space (lower 
graph)R. 



^ Note that in the abscissa axis of the Fourier transform we 
indicate the period rather than the frequency 
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This signal oscillation represents an estimate of 
the thermally-induced signal instability caused by the 
Planck Sorption Cooler working with degraded perfor- 
mances. It is worth underlying that this estimate was de- 
rived from simulations^ which do not represent a quan- 
titative prediction of the expected Sorption Cooler per- 
formances, but were aimed at understanding qualitatively 
the impact of compressor non-idealities on the tempera- 
ture stability. This oscillation pattern is used in our paper 
merely as an example of the application of our analysis to 
signal oscillations of complex shape. 

From the figure it is apparent that the fluctuation is 
dominated by two main periods (667 s and 4000 s) but 
it also contains a high number of harmonics with periods 
down to about 40 s. The inset in the lower graph of Fig. ^ 
is a close-up of the high frequency tail of the fluctuation 
spectrum. 



By applying Eq. ( pij[ ) to the Fourier transform of the 
signal we are able to calculate the expected peak to peak 
amplitude of this effect when it is processed by our de- 
striping code and projected onto a 13.7' map. In the case 
discussed here the computation yields an expected value 
of 13.37 fiK for the flnal peak-to-peak amplitude, which 
is in excellent agreement with the value obtained by the 
standard mapping procedure (see Fig. 

Note that in Fig. || the morphology of the residual 
systematic on the map (which is relative to the signal os- 
cillation shown in Fig. ^ is very similar to the morphology 
shown in the bottom map of Fig. |4| (relative to a "pure" 
periodic signal with T{ = 4000 s) . This is the consequence 
of the fact that the oscillation in Fig. ^ is dominated by 
the 4000 s wave (as shown in the Fourier transform) which 
determines the main visible structure in the map. 
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Fig. 8. Map (pixel size = 13.7') after destriping of the 
signal fluctuation plotted in Fig. 0. 
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Fig. 7. Signal oscillation in time domain (upper graph) 
and Fourier space (lower graph) representing a simula- 
tion of the thermally-induced signal instability caused by 
a Sorption Cooler working with degraded performances. 
The inset in the lower graph is a close-up of the high fre- 
quency tail of the fluctuation spectrum. 



^ Simulations were run using the Sorption Cooler thermal 
model implemented with the SINDA simulation software 



The same procedure was also applied to other two sig- 
nal fluctuations with very different shapes, which repre- 
sent thermally-induced fluctuations for various Sorption 
Cooler performance scenarios. The results of the three 
tests summarised in Table |l| show a very good agreement 
between the final peak-to-peak amplitudes obtained di- 
rectly from maps and calculated with Eq. ([lo| ) , which con- 
firms the general validity of our approach. 

The values reported in the last two columns of Table |^ 
are relative to 30 GHz LFI maps with a pixel size 0pixoi = 
13.7'. Because the optical beam size will be approximately 
three times larger than 0pixoi, the actual amplitude of the 
systematic effect on the instrument angular scale would 
be much smaller, ranging from ~ 0.5 /uK to ~ 5 /iK for 
the three cases reported in Table 0. 

Although further analysis (which is out of the scope 
of our paper) is needed to understand the impact on the 
recovered amplitude and polarisation power spectra, our 
analysis can give rather precise estimates of the damping 
factor of the peak-to-peak amplitude of the fluctuation 
at the instrument output after the projection in the flnal 
maps. 
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Peak-to-peal< amplitude on 
final maps after destriping 
with 13.7' pixel size 


Shape of 

signal 
fluctuation 


Sorption 
Cooler 
performance 


Peak-to- 

peal< 
amplitude 

(mK) 


From 
analytical 
calculation 


From 
simulation 




Nominal 


1.06 


2.84 


2.31 




Degraded 


4.30 


1337 


13.04 




Highly 
degraded 


8.24 


22.88 


21.65 
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Table 1. Application of our analytical procedure to sig- 
nal fluctuations of different shape and amplitude. These 
fluctuations represent the effect of the Sorption Cooler 
working with different levels of performance. 



4. Conclusions 

In this paper we have proposed an analytic method to 
estimate the effect of periodic signal oscillations on CMB 
maps to be produced by sky scanning experiments. When 
these signal fluctuations are processed with a destriping 
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